PARTIAL LEGENDRE TRANSFORMS OF NON-LINEAR EQUATIONS 



PENGFEI GUAN AND D. H. PHONG 

Abstract. The partial Legendre transform of a non-linear elliptic differential equation 
is shown to be another non-linear elliptic differential equation. In particular, the partial 
Legendre transform of the Monge- Ampere equation is another equation of Monge- Ampere 
type. In 1+1 dimensions, this can be applied to obtain uniform estimates to all orders for 
the degenerate Monge- Ampere equation with boundary data satisfying a strict convexity 
condition. 



1. Introduction 

The maximum rank and regularity properties of degenerate fully non-linear equations are 
still largely unexplored, despite their considerable interest for many geometric problems. 
For example, it is still an unresolved problem raised by Donaldson [9J [10] to determine the 
precise regularity of geodesies in the spaces of Kahler potentials and of volume forms on a 
Riemannian manifold. These are given respectively by solutions of a degenerate complex 
Monge- Ampere equation and an equation introduced by Donaldson [10J. Many existence 
and regularity properties have now been established for these equations (see e.g. [H HU [5j 
[151 EH El US HH] and references therein), but it is not known how close they are to optimal. 

In |14j . the maximum rank property has been established for several special cases of the 
degenerate real Monge- Ampere and Donaldson equations, for Dirichlet data satisfying strict 
convexity conditions. Thus the natural question arises of whether the regularity of solutions 
of fully non-linear elliptic equations can be established, assuming that it is already known 
that they have maximum rank. 

A potentially useful feature of the maximum rank property is that it allows the use of a 
partial Legendre transform. In fact, the partial Legendre transform was already exploited 
by D. Guan [12j in determining geodesies for the space of Kahler potentials on toric varieties, 
and by P. Guan [13] and Rios, Sawyer, and Wheeden [201 El] in their study of the local 
regularity of certain degenerate Monge- Ampere equations. A first goal of this paper is to 
refine their analysis and show that, even though the partial Legendre transform / of a 
function u is not a local expression of u, the partial Legendre transform of an elliptic PDE 
in u is another PDE in / which is again elliptic. In particular, the original Monge- Ampere 
equation can be transformed globally into another dual equation, again of Monge- Ampere 
type, but which does not seem to have been encountered before in the literature and may 
be of independent interest (see Theorem [H (b)). The partial strict convexity properties of 
one equation are then equivalent to C 2 estimates for its dual, and one can expect a more 
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general correspondence between bounds for their derivatives. The second goal of this paper 
is to apply this principle in the simplest case of the 1 + 1 real Monge- Ampere equation 
(which coincides with the 1 + 1 Donaldson equation). As a consequence, we obtain C°° 
bounds for this equation which depend only on the Dirichlet data, and in particular which 
remain uniform as the equation degenerates (c.f. Theorem [2]). We note that in this case, 
by |12| . the solution of the limiting equation is already known to be smooth, so the real 
interest of the result lies in the uniform validity of the approximation. 

2. Legendre transforms 

The main goal of this section is to work out the partial Legendre transforms of fully 
non-linear elliptic PDE's in some generality. We shall find that, just as in the case of the 
full Legendre transform, they are given by elliptic PDE's. 

2.1. The full Legendre transform. We begin by re-visiting the standard Legendre trans- 
form. Let u(x) be a strictly convex function on R". Then u(x) defines a Legendre change 
of variables 

du . . 

(2.1) X ^y=— {x ). 

Clearly, the Jacobian of this change of variables is = ^.q x = Ujk- The strict convexity 
implies that the map x — > y from R n to its image is invertible. The Jacobian of the inverse 
y — > x is given by the inverse u^ k of the Jacobian Ujk of x — > y. Associated to the function 
u(x) is also its Legendre transform f(y), defined by 

du 

(2.2) f(y)=xy-u(x), y = — . 

Differentiating this relation with respect to y shows that the Legendre change of variables 
defined by f(y) is the inverse map y — > x 

(2.3) y -+ x = §£, 

dy 

and we have the following exact analogues of the earlier formulas for u, = Q x .Q Xk = u^ k . 
A partial differential equation of the form Ffuj^) = can be viewed as a partial differential 
equation in /. Its linearization has principal symbol = 'f kq £,p£,q- Thus the ellipticity 
of the equation in u implies the ellipticity of the equation in /. In particular, a Monge- 
Ampere equation for u, 

(2.4) det( 7 — — ) = K 

OXjOXk 

is equivalent to a Monge- Ampere equation for / 

(2.5) det{ ^f ) = K -\ 

We note that the changes of variables x — > y and y — > x in (|2.1l I2.3() are unaffected if 
u and/or / are shifted by independent constants. The relation (|2.2p can be viewed as a 
canonical way of fixing the relative normalization of / and u. 
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2.2. The partial Legendre transform. We come now to the situation of main interest 
to us, namely partial Legendre transforms of functions u(x, t) which are periodic and satisfy 
a a strict convexity condition in x. More specifically, let be the basis vectors for R n , that 
is, ej has component 1 in the i-th position, and all its other components are 0. We consider 
functions u(x, t) on R" X I, I = (0, 1), satisfying the periodicity condition 

(2.6) u(x + ei) = u(x), x G R n , l^i^n, 

and the strict convexity condition 

(2 - 7) d^k +6 ^ >0 - 

Thus u can also be viewed as a function on X x 7, where X is the torus X = (R/Z) n . 
We define the following Legendre change of variables 

du 

(2.8) x^y = (y fe ), y k = — + Xk . 

The inverse map y — >■ x is well-defined and unique as a map from R n to R n . To see this, 
set v = u + ^|x| 2 , and note that v — > oo as x — > oo. Thus v admits a minimum. After a 
translation if necessary, we may assume that this minimum is at 0, and in particular, that 
Vi> vanishes there. It suffices then to see that, for any y, the equation 

dv 

(2.9) —( Xt ) = t y k , O^t^l 

can be solved by the method of continuity: the Jacobian (vjk) is invertible implying openness 
of the set of t's for which the equation is solvable. And since \Vv\ — > oo as \x\ — > oo, we 
also find that \x t \ is bounded, and the set of such t's is also closed. 
The maps x — > y and y — > x satisfy the following transformation laws, 

(2.10) y k (x + ei) = y k (x) + S ik , x k (y + a) = x k (y) + S ik . 

The condition for y k follows immediately from its definition and the fact that u(x) is peri- 
odic. To establish the condition for y, just observe that 

,„ _. . . du du 

(2.11) y k + d ik = q^( x ) + %k + o ik = Q^( X + e i) + i x + e i)k 

and the assertion follows by the uniqueness of the inverse map y —> x. 
Clearly, the Jacobian of the map x — )■ y is given by 

(2 12) — - d2u |5 =„.. 

dxj dx k dxj jfc 

Consequently, the Jacobian of the inverse map y — > x is given by 
(2.13) 

oyk 

All these expressions are periodic, and descend to equations on the torus X. 
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So far, we have discussed the Legendre maps defined by the function u. We now define 
the Legendre transform / of u itself by the following formula, 

j Qix 

( 2 - 14 ) f(y) = —j\x - y\ 2 - u(x), y i = ~fa: + x i 

for y 6 R n . We note that, in view of the transformation laws fj2. lOf) for x and y under period 
shifts, the function f(y) is actually periodic, and thus can be identified with a function on 
the torus X. Again, the inverse map y — > x can be viewed as the map associated with the 
function /, 

df 

(2.15) y ' ^ ■ >l, 

and its Jacobian is given by 

(2-16) P- = J-jt- + S jk = g> k = h jk . 

oyj oyjdy k 

In particular, the Legendre transform / satisfies the strict convexity condition 
(2.17) D 2 y f + I>0 

and we still have the relation 

Consider now the change in variables R n x I — >■ R n x I defined by 

du 

(2.19) (x,t) ->• (y,s), yj = — + Xj, s = t. 

The Jacobian of the inverse change of variables is given by 

(2.20) t s = l, t yp = 0, (x k ) yp =g kp , (x k ) s = -g kp u Xpt . 

It follows that the rule for differentiating a function F(x, t) with respect to the variables 
(y, s) is given by 

(2.21) F Vp = F Xj J', F s = -F Xj g^u Xpt + F t . 

The dependence of the partial Legendre transform on the additional variables s and t is 
now conveniently described by the following three equations 

(2.22) dy.ut = -d s Xj, d s u t = K (detg)~\ f s = -u t . 
Here K is the (n + 1) x (n + 1) determinant det(D^. t u + I x ), 

(2.23) K = det (D 2 xt u + /,) = u tt (detg) - G^u tX] u tXp , 

where G 3P = (detg)g JP is the matrix of co- factors of the metric gij. To see the first equation, 
we compute both sides. On the left, we have d Vp ut = ut Xj {xj) Vp = ut Xj gi p . On the right, 
we have —d s x p = —(—g P3 u tXj ) = g p ^u tXj also, as required. Next, we apply the rule for 
differentiation of the previous paragraph and obtain 

(2.24) d s u t = -u tx .g> p u tXp + u tt = (detg)~ l (u tt detg - G ]P u tX] u tXp ) 



as claimed. Finally, differentiating the defining formula (|2.14j) for / gives 

(2.25) d s f = (y - x) ■ x s - (u x ■ x s + u t ) = (y - x - u x ) -x s -u t = -u t . 

All three identities in (|2.22p have been proved. They imply readily the following two iden- 
tities, which we also need later 

(2.26) u tXj = -f syk h kj , u tt = -f ss + f syj fsy k h jk . 

2.3. The partial Legendre transform of non-linear PDE's. We consider now a fully 
non-linear equation of the form 

(2.27) F{D 2 u) = 

on Ix/, where the unknown u is required to satisfy the strict convexity condition D 2 u+L x > 
0, and the equation is assumed to be elliptic. We would like to view this equation as an 
equation for the Legendre transform / of u. Note that / is a non-local quantity in u. 
Nevertheless, we have 

Theorem 1. Let X = (R/Z) n be the n- dimensional torus, and let I = (0,1). Let u(x,t) 
be a function on X x L satisfying the strict convexity condition \2. 7\ ). Let (x, t) — > (s = t, y) 
be the partial Legendre transform as defined by A2.3\) . and let f be the partial Legendre 
transform of the function u as defined by \2.1J$ . 

(a) If F(D 2 u) = is a second-order elliptic PDE for u, then it can also be viewed as a 
second-order elliptic PDE for the partial Legendre transform / of u. 

(b) If particular, u satisfies the Monge- Ampere equation 

(2.28) det (D 2 x t u + L x ) = K 

if and only if its partial Legendre transform / satisfies the following equation onlxl also 
of Monge- Ampere type, 

d 2 f 

(2.29) J +Kdet ( D 2 f + I ) = . 



Proof. From the discussion in the preceding section, the partial Legendre transform / of 
u is a well-defined function on X x /, djdkU + 5jk are given by the inverse of the matrix 
djdkf + 6jk, and ut Xj and uu are given by the expressions (|2.26[) in the second derivatives 
of /. Thus the equation F(D 2 u) is automatically a second order non-linear equation for /. 
To verify the ellipticity of the equation viewed as an equation for /, we work out first the 
linearized operator of F(D 2 u), keeping variations in Su, 

dF dF dF 

(2.30) 5F = -T—(5u) tt + (5u) tX] + {5u) XjXk . 

du tt du tXj 3 du XjXk 3 

We need to express this quantity in terms of the derivatives of Sf. In view of the expression 
(|2.26p for Suft and 5ut Xj , we have 

Su tt = -5f ss + 25f sy J syk W k - f sy J syk 5fi k 

(2.31) 5u tXj = -5f syk h k i +f syj 5f k , 
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where the indices are raised or lowered using the metric hji~. Thus we have 



(2 ' 32) + <-s£w~+£;/«- 5 ^w'*. 

This means that, if r and £j are respectively the variables dual to t and Xj, the symbol 
cr(r, £) of the linearized operator is given by 

duu du tXj du tt 

dF dF dF 

(2-33) +(jr—fsy j fsy k - /sj/* + q- )£i£fc- 

We shall show that <t(t, £) is positive definite if the original equation F(D 2 u) = is elliptic. 
Introduce the variable m = hP k ^ for convenience. Then completing the square in r gives 



dF 1 dF dF Vj 



du tt 2du tXj du tt yj [W 

V du tt 

/ ^ 1 r, dF dF , dF 

To prove part (a) of the Theorem, it suffices then to show that that the second term on 
the right is strictly positive for r] ^ when F is elliptic. Since the symbol of the linearized 
operator when the unknown is u is given by 

dF 2 OF ^ dF 

(2.35) — - r 2 + rij + — £& k 

du tt dutxj du x . Xk 

its ellipticity does imply the positivity of the second term on the right hand side of ()2.34p . 
Part (a) is proved. 

Part (b) follows immediately from the identities in (|2.22|) . 

(2.36) f„ = -d s u t = -Ki&etg)- 1 = -Kdet (D 2 y f + /). 
The proof of Theorem [T] is complete. Q.E.D. 

We conclude this section with a few remarks. 

• Part (b) of Theorem [1] can be viewed as a refinement of several earlier results in the 
literature using the partial Legendre transform: when K = 0, it reproduces the result of 
D. Guan |12j . For general K, and when the considerations are local (instead of on a torus 
as here), then P. Guan [13] in dimension n = 1 and Rios- Sawyer- Wheeden |20j for general 
dimension n have shown that the coordinates Xj, viewed as functions of (y,s), satisfy the 
following elliptic, non-linear system of equations 

(2.37) d 2 Xj + JL {Kdet{ pL )) = , 
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This system of equations follows immediately from differentiation of the equation (|2,29p 
with respect to yj. 

• The presence of the background symmetric form 5jk is very similar to the presence of 
the Kahler form for the complex Monge- Ampere equation (0 + ^dd(p) n = F0 n . 

• The correspondence between an equation in it and its "dual" equation in / can provide 
non-trivial information. For example, it is not evident that the equation (|2.29p admits 
smooth solutions for given Dirichlet data, even when K is a strictly positive constant. On 
the other hand, the existence of such solutions is an immediate consequence of the existence 
of smooth solutions for the dual equation (|2.28p . which can be established by the theory of 
Caffarelli-Nirenberg-Spruck [8], with the improved barrier arguments of B. Guan 

Similarly, lower bounds for D 2 u + / are equivalent to upper bounds for D 2 f + / and vice 
versa, and the problems of partial C 2 estimates and partial strict convexity are in this sense 
"dual" . For example, the C 2 estimates for the original Monge- Ampere equation (I2.28h can 
be established by traditional methods as in [8], or as a consequence of the convexity results 
for the dual equation (12.29p . using for example the recent results of Bian-Guan [U [2]. 

• Beyond the Monge- Ampere equation, the partial Legendre transforms of Hessian equa- 
tions may be of interest. By (|2.26p . 



-u syi X 1 1 A x x 



f-K(detg) 1 -u^A^ ■■■ -u^A^ 

(2.38) Dl it u + 1 = 

V -u syn \- 1 A" 1 J 

where Aj are the eigenvalues of g = (uij + 5ij), and we have denoted for convenience all 
quantities associated with the partial Legendre transform by a ~ (e.g. / is now denoted u, 
and K is the Monge- Ampere determinant of D 2 s u + /). The standard formula for the k-th 
symmetric function of the eigenvalues of a matrix is 



k\ 

We apply this formula to the the symmetric function a of the n + 1-dimensional matrix 
Dl t u + I. We find 



k-l 

(2A0)a k (D 2 Xjt u + J) = akir 1 ) - £ Ar 1 ■ ■ ■ u^.Ar 1 + ^(det^)- 1 )- 

h—ik-x ij=l 
From here, it is easily seen that the Laplace equation uu + An = K\ gets transformed into 
(2.41) K + K l det~g = a n _ 1 {g) 

For general k, the equation ak(u) = gets transformed into 

k-l 



We note that when k = (n + 1), the above identity recovers (|2.29l) . 
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3. The 1 + 1 Monge- Ampere equation 



In this section, we consider more specifically the case n = 1 of the Monge- Ampere equa- 
tion. In this case, the equation becomes the following equation on X x I, 

(3.1) u tt (l + u xx ) - u 2 xt = e 

where X = R/Z is a circle and e > is a constant. The solution u is required to satisfy 
D 2 t u + I x ^ 0, and we impose the Dirichlet condition u(t,0) = u°(x), u(t, 1) = u (x), with 
u l G C°°(X) satisfying the strict convexity condition, 

(3.2) u xx + 1 > A 
for i = 0, 1 and A > a constant. We have 

Theorem 2. Let u be the solution of the equation Ii3. 1)) on X 1 x I, with D xt u + I x Si and 
smooth Dirichlet data satisfying A3.2\) . Then for any non-negative integer N , there exists 
M(N) and a constant CV depending only on A > and the C M ( N ^ norms of the Dirichlet 
data u° , u 1 so that 

(3.3) Y. H^M| c o (x i xJ) < C N . 

a+b^N 

In particular, the constants Cm are independent of e. 

The rest of this section is devoted to the proof of Theorem [2j To apply the partial 
Legendre transform as in §2, we need the strict partial convexity of u. This follows from 
Theorem 1 of |14| . But since the present case is particularly simple, we can supply the short 
proof for the convenience of the reader: 

Lemma 1. Let u(x,t) be the solution of the equation \3. 1}) . as specified in the statement of 
Theorem^ Then u xx (x,t) + 1 ^ A for all x G X. 

Proof of Lemma\^ Let Ao = min^ j)eXxi( u xx + 1), and set if(x,t) = u xx + 1 — Ao- We 
establish a strong maximum principle for ip. If (p vanishes on the boundary, the lemma 
is proved. We shall show that in a neighborhood of any interior zero of ip, ip satisfies an 
elliptic differential inequality equation of the form 

(3.4) F l ^ l3 ^C\Vip\. 

Here we have denoted by F the function F{D^. t u + I x ) with F(M) = det(Mjj), for M any 
symmetric 2x2 matrix, and i™ = gj^.. - The constant C is required to be independent 
of the point (x,t) but may depend on everything else. By the strong maximum principle, 
this would imply that <p vanishes in a neighborhood of any interior zero. If the set of such 
zeroes is not empty, then tp vanishes identically. By continuity ip would again vanish on the 
boundary, and the lemma is proved in all cases. 

The equation (|3,ip can be written as F(D^. t u + I x ) = e. Differentiating the equation 
successively gives 



(3.5) F ' Uij x — 0, F ^Uij xx + Uij x U};\ x — 0. 
Thus F^ipij = —F l i ,kl Uij x Uki x: and more explicitly, 

(3.6) F lj ipij = -2u ttx u xxx - 2u 2 xxt = -2u ttx ip x - 2tp x . 
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The inequality {33} follows. Q.E.D. 



Let f(y, s) be now the partial Legendre transform of u, as defined in section. When 
n = 1, the equation for / simplifies to 

(3.7) Lf = f ss +ef yy = on X x I, 

and f(y, 0) and f(y, 1) are given by the Legendre tranforms f°(y) and f 1 (y) of the functions 
u°(x) and General linear elliptic theory says that any derivative of / can be bounded 

in terms of the boundary data and the ellipticity constant e. However, we require estimates 
which are uniform as e — > 0, and such estimates do not seem to have been written down in 
the literature. We provide below a brief and explicit derivation of estimates uniform in s, 
exploiting the simple form of the equation and of the boundary in the present case. More 
precisely, we shall establish the following lemma: 



Lemma 2. Consider the Dirichlet problem for the Laplacian L in (3.1) with e a constant 
satisfying < e < 1. Then for any m, k we have 

where C m ^ are constants which depend only on the Dirichlet data (and on m and k). In 
particular, C m ^ are independent of e. 

Proof. Clearly d^d^f satisfies the same Laplace equation. By the maximum principle, we 
have then 

( 3 - 9 ) Hd™ds/llc°(Xx/) = W d y' d sf\\c°(XxdI)- 

We shall show the right hand sides can be estimated in terms of the Dirichlet data alone, 
for arbitrary m and k = or k = 1. Assuming this, it follows that the left hand side is 
also bounded in terms of the Dirichlet data alone for these values of m and k. Since the 
equation implies that d^d^f = — £<9™ +2 <9g ~ 2 f for k ^ 2, it follows that uniform bounds for 
dj^dgf for arbitrary k follow from the special cases k = and k = 1, and the lemma would 
be proved. 

We return to the proof of bounds for ||<9™c^ f\\c°(Xxdl) wnen A; = or fc = 1. When 
k = 0, they are obvious, so we concentrate on the case k = 1. Let / be a function with the 
same boundary values as / (e.g., / = tf(y, 1) + (1 — t)f(y,0)). Set 

(3.10) w(y,s) = -As 2 + Bs 

for constants A, B > which we shall choose in a moment. Since L(d y n f) = 0, we have 
(3.H) L(d?U-f))-L{8yf)><* 

where c\ is a constant depending only on the Dirichlet data and independent of < e < 1. 
On the other hand 

(3.12) Lw = -2A 

so we can choose A large enough so that Lw ^ L{d y n {f — /)) on X x I. Now d y n (f — 
f) vanishes identically on the boundary, so if we choose B large enough so that w ^ 
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everywhere, we shall also have d™(f - /) ( id on the boundary. Thus, by the comparison 
principle, <9™(/ — /) ^ w on X x I. Since both of these functions vanish at s = 0, we obtain 

(3.13) 8 s d™f < -2As + B + d s 0™f, 

which is an upper bound for d s d!^f. Applying the argument to — / instead of /, we obtain a 
lower bound for d s d™f at the boundary points (y, 0). Since the argument at the boundary 
points (y, 1) is identical, the proof of the lemma is complete. 

Next, we show that C N bounds for x (viewed as a function of (y,s), together with a 
strict partial convexity bound, imply C M bounds for the original function u: 

Lemma 3. For any non-negative integer M , we have 

(3.14) \\d™d b t u\\ CO{x t xI) ^ C M 

m+b^N 

where Cm is a constant depending only on the Dirichlet data it , it 1 , the lower bound A > ; 
and the C° norm of a finite number N(M) of spatial derivatives d™x of the function x. 

Proof of Lemma [3J First, we show that bounds for d™x imply bounds for d x u. This is 
an easy consequence of the following formula, which is itself a consequence of the chain rule 
established in the previous section: 

y (u xx + i) m (u xx + iy m 1 

where P is a generic notation for a polynomial in all its entries for all m ^ 2. 

Now for m = 2, the bounds of d^u in terms of boundary data alone are a special case 
of the C 2 estimates for the Dirichlet problem for the Monge- Ampere equation [8]. Note 
that these bounds do not require a strictly positive lower bound for the Monge-Ampere 
determinant, and thus give bounds which are uniform in e in our case. 

Assume that bounds depending only on the Dirichlet data and a strictly positive lower 
bound A for u xx + 1 have been established for m. In view of the above formula for d™x, 
it follows that such bounds for reduce to such bounds for d™x on X x I. By the 

maximum principle for d™x, this reduces in turn to bounds for <9™x only on the boundary. 
But the same formula (|3.15p for higher derivatives above shows that d™x on the boundary 
X x dl is determined completely by the boundary data. This establishes the desired bounds 
for ||<9>|| c o (X x 7) . 

Next, we consider mixed derivatives of the form d x dtu. It is again easy to establish the 
following general formula linking d™d s x and d x dtu for all m ^ 0, 

V ' y (n^. + l) m+1 (u xx + l) 2m +! 

where P is again a polynomial in all its entries. In view of Lemma [21 the left hand side can 
be bounded uniformly in terms of the Dirichlet data and the lower bound A for u xx + 1. 
Thus the formula implies that dJ^^dtu can be similarly bounded if dJ^dtu is. Since u x t 
is bounded by the Dirichlet data in view of the C 2 estimates, we obtain by induction the 
uniform boundedness of d^dtu for all m. 
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Finally, by differentiating the original Monge- Ampere equation, we can show inductively 
on the number 6 of i derivatives in d x d b u that they are in turn bounded. First, we show 
this for 6 = 2. Differentiating the equation m times with respect to x gives 

(3.17) d™u tt = l — P{d 2 x u, • • • , d™ +2 u, d t d x u, • • • , d t d™ +1 u, 8 2 t u, • • • , d 2 t dZ l ~ l u) 

with P a polynomial in all its entries. Since d™uu is bounded for m = by the C 2 estimates, 
and d%u, d x dtu are now known to be bounded, the formula allows us to show by induction 
on m that d™dfu is bounded. Next, assume that d^d^u is bounded for an integer 6^2 
and all non-negative integers m. We shall show that this remains true if 6 is replaced by 
6+1. Differentiating the equation 6—1 times gives 

(3.18) d b t +1 u = — l —p(d x d t , • • • , d x d b u, d 2 x d t u, ■■■ , d 2 x d b t u) 

"ix t l 

where P is a polynomial in all its entries. This shows that d^d^^u is bounded for m = 0. 
Thus it suffices to show that if d x d b+l u is bounded for all non-negative integers £ ^ m, then 
the same is true if m is replaced by m + 1. But this follows readily by differentiating the 
previous formula, 

(3.19) d ™+i d b+i u = d ™+i{_L_ P(0 x9t ,... ,d x d b t u,d 2 x d t u,--- ,8 2 x d b u)}. 

Since the right hand side involves terms with at most 6 derivatives in t, and all derivatives 
in x of such expressions are bounded, the desired bound for follows. Q.E.D. 

Putting together Lemmas 2-4, we obtain Theorem [2j 
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